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Abstract 

We study a similarity transformation to construct an effective Hamiltonian sys- 
tematically, which does not contain particle-number-changing interactions, by means 
of Fukuda-Sawada-Taketani-Okubo's method. We show that such Hamiltonian can be 
constructed from Feynman diagrams and give rules for constructing it in the Light- 
Front Yukawa model. We prove that it is renormalized by the familiar covariant per- 
turbative renormalization procedure. It is very advantageous that the effective Hamil- 
tonian can be obtained from our rules immediately. We also numerically diagonalize it 
to the second order in the coupling constant as an exercise. 



* E-mail: yukilscp@inbox.iic.kyushu-u.ac.jp 



1 



typeset using pTpTEX.sty <ver.l.O> 



1. Introduction 



For many years, the bound state problem for QCD has been studied very well but has 
not come to the complete understanding until today. The main difficulty is that the low 
energy djTiamics of QCD needs nonperturbative calculations. One of the convenient tools 
to solve the problem is a Hamiltonian approach in Light-Front (LF) field theory which is 
quantized on the equal LF time surface. 

An important feature of the LF field theory is that the vacuum is trivial; the Fock 
vacuum of the free part of the Hamiltonian is the true one. It is very useful in the relativistic 
bound state problem because we do not need to worry about solving the complex vacuum 
in contrast to the usual equal-time (ET) field theory. Even though the vacuum is trivial, to 
solve the Schrodinger equation for the relativistic bound state is difficult because it is natural 
that the eigenstates are constructed by the superposition of an infinite number of particle 
states allowed by symmetries of the Hamiltonian. Tamm-Dancoff (TD) approximation. Ill 
which truncates the Fock space, that is, limits the number of particles concerned with the 
interaction, simplifies the practical calculations. While TD approximation was originally 
proposed in the ET field theory. Perry, Harindranath, and Wilson suggested a LF version of 
it.& They pointed out that since the truncated states do not consist a complete set of the 
Hamiltonian, ultraviolet (UV) divergences are nonlocal and noncovariant, and counterterms 
which should cancel the divergences depend on the sector of the Fock space within which 
they act. This is called the problem of "sector-dependent counterterms". 

A similarity transformation of the Hamiltonian does not change the eigenvalues and is 
useful for getting the effective Hamiltonian. There are two types of it. One is the transfor- 
mation in momentum space, which is equivalent to integrating out states which exchange 
energies more than some energy cutoff, proposed by Glazek and Wilson, & and indepen- 
dently by Wegner.& It is interesting that it gives the nonperturbative low energy physics 
and a logarithmic confining potential in LF QCD,i) although it is hard to get the effective 
Hamiltonian even in the lowest order in the coupling constant. The details of this method 
and applications are discussed in Ref. and recent progress is seen in Ref. |^). 

Another is the transformation in the particle number space, which reduces the Hamilto- 
nian to one which has no particle-number-changing interactions so that the transformed one 
can be solved easily, and was proposed by Harada and Okazaki.i' It can avoid the problem of 
sector-dependent counterterms because the origin of it is the general property of the Hamil- 
tonian that it has particle-number-changing interactions and its eigenstates need an infinite 
number of particles. But its actual calculations are complicated and tedious. Although their 
method was considered in the LF field theory, it is not new in the ET context and has been 
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used for getting the TMO potential of nuclei, u We call it Fukuda-Sawada-Taketani-Okubo's 
(FSTO's) method. ©'0' This method gives us an easier way for constructing the effective 
Hamiltonian and seems to be promising in the LF framework. However, it lacks manifest 
Lorentz covariance, and therefore it is difficult to tell what sort of divergences the effective 
Hamiltonian has before doing actual calculations. It makes the renormalization procedure 
more complicated than the usual covariant perturbation theory. It is highly desired to make 
transparent how the divergences emerge in the FSTO's framework. 

The purpose of this paper is to show that the effective Hamiltonian constructed by the 
FSTO's method can be immediately obtained from Feynman diagrams as the S-matrix ele- 
ment can be in the covariant perturbation theory, and that one can use the usual renormal- 
ization procedure to renormalize it. It makes the construction of the effective Hamiltonian 
systematic and easier. Especially, it allows us to perform higher order calculations. 0) 

Our strategy is the following. First, we show that the FSTO's effective Hamiltonian is 
a sum of the auxiliary operator G and its products to the fourth order of the interaction 
part of the Hamiltonian. The advantages of using G are that it is constructed by the same 
Feynman diagrams as those for S-matrix elements, and that it has no particle-number- 
changing interaction so that one can easily calculate normal-ordering of the products of 
those. We give the rules for constructing G from Feynman diagrams. They are a little 
bit differ from the familiar Feynman rules in the covariant perturbation theory. The set of 
those rules is one of the main results of the present paper. Our method is more convenient 
and powerful for constructing the effective Hamiltonian and discussing the renormalization 
of it than the other similarity methods. Then, we find that there are three types of UV 
divergence. One is the familiar loop divergence which can be renormalized by the usual 
renormalization procedure. The second comes from the difference between our construction 
rules and Feynman ones. We can make it harmless by using an ambiguity. The above two 
only emerge in G. The last emerges from the products of the renormalized G"s. We show 
that it works as a box counterterm, which is needed to cancel the cutoff dependence of 
the eigenvalue, in diagonalizing the effective Hamiltonian. 

This paper is organized as follows. In Sec. |], we briefly review the FSTO's method. 
In Sec. 0, we show that the effective Hamiltonian constructed from the FSTO's method is 
written in terms of the auxiliary operators F or G. F is written in terms of G. We consider 
the LF Yukawa model, in which Lagrangian and Hamiltonian are given in Appendix ^ as 
a concrete example, and give the rules for the construction of G from Feynman diagrams and 
how to renormalize it and the effective Hamiltonian. In Sec. ^, we summarize and discuss 
the validity of our method. In Appendix as an exercise, we calculate the eigenvalue of the 
ground states of the effective Hamiltonian up to the second order in the coupling constant in 
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the LF Yukawa model. In Appendix |C[ we explain an ambiguity of the energy integrations. 



§2. Review of the FSTO's method 

This section briefly reviews the FSTO's method partly following Ref. |ll]). The FSTO's 
method is to reduce a Hamiltonian to the block-diagonal form using a similarity transfor- 
mation. We can get the effective Hamiltonian for the subspace of the Fock space without a 
loss of the necessary information. 

We want to solve the Schrodinger equation 

H\^)=E\^), (2-1) 

for the second-quantized Hamiltonian H which consists of the free part Hq and the interaction 
part H': 

H = Ho + H'. (2-2) 

The eigenstate jip') can be expanded by the complete set of Hq. To divide the Fock space 
into two, we introduce a projection operator r] which commutes with Hq: 

V\^) = iV'i), .2 3) 

{1-vm = 1^2), 

where both {ipi) and \1IJ2) are the states written in terms of the eigenstates of i^o- iV^i) is 
the state in the target Fock space at our disposal. For our purpose, we restrict t] to the one 
which selects the states with definite number of particles. In the matrix notation, rj and \^) 
are written as 

'=(::)■"•=(::::)■ 

respectively, and we express an arbitrary operator O as 

Q^(On Ou\ Ou = vOv, 0^2 = vOil-v), .2-5) 

[O21 O22)' O21 = (1 - r/)Or/, 022 = {l-r])0{l-ri). 

Let us consider the similarity transformation that leads the Hamiltonian to the block- 
diagonal form in the following way 



1^2)' / I * M 1^2)' 



(2-6) 

where "*" abbreviates the part which is not needed for our purpose, and 

(2-7) 



1^1)' \ 1^1) 
1^2)' j I 1^2) 



4 



\ipi)' completely decouples from \ip2)' ■, so that we can concentrate only on the equation for 
ifeff and I'j/'i)' in the subspace selected by rj: 

i/effl^l)' = ^l^l)'. (2-8) 

We set an ansatz for the similarity (unitary) transformation operator 

\U2i U22 I [A{l + AU)-y^ (l + AAt)-i/2 ]■ ^ ) 

It is convenient to use J 

instead of A. In order for the off-diagonal parts in Eq. ( p-6| ) to become zero, J must satisfy 

(1 - r]) {H'J + [Ho, J] - J{H'J)) r] = 0, (2-11) 
where we have introduced the notation ( ) as 

{0)^vOv + il-v)Oil-v)=i ° V (2-12) 

y (J U22 J 

Using Eqs. (|2-6| ), (|2-9| ), and ( |2-11| ), we formally obtain the effective Hamiltonian 



H,s = Ti{J^J)-''\J^HJ){j'^J)-^/\ (2-13) 

but J is unknown so far. 

In the interaction picture, Eq. ( P-11 ) may be regarded as a differential equation 



(1 - ^)^^^ = (1 - ^) \H\t)J{t) - J{t){H\t)jm V, (2-14) 

where J(0) is identified with J. H'(t) and J{t) are defined in this picture as 

H'{t) = e'"'''H'e-'""'e'\ J{t) = e'^"' J e''""' e'\ (2-15) 

where e*^* is an adiabatic factor with e 0+. We hereafter omit it in order to keep equations 
simple. Provided that we set the initial condition 

J(-oo) = 1, (2-16) 

J can be found order by order. Alternatively, if we solve 

= H'{t)V{t) - V{t){H'{t)V{t)), (2-17) 



under the initial condition 

V{-oo) = 1, (2-18) 

J is given by 

It is easily found that the solution of Eq. ( |2-17| ) is given by 

V{t) = U{t){U{t))-' 

U2l{t) U22{t) ) [ 

U2imii\t) 1 




(2-20) 



where U{t) is the usual time evolution operator 

U{t) = Texp |-i y* dt'H'it')^ . (2-21) 

Uiiit) and U22^{t) are the inverse operators of Unit) and U22(t) respectively. J is explicitly 
obtained from Eqs. ( |2-19D and (2-20), and then we expand it in the order of H' to obtain 
the effective Hamiltonian perturbatively: 



oo 



J = I ^ , ° I = 1 + y (2-22) 

where n is the order of H'. The resultant effective Hamiltonian becomes 
iJeff = Iv + H' + H'Ji + H'J2 

+H'J3 + ^[jUi, H'Ji]^ T] + h.c. + 0{H"), (2-23) 

to the fourth order. It is known that this can produce the TMO potential in the symmetrical 
pseudoscalar pion theory with the pseudovector coupling, which explains the properties of 
the deuteron very well. 

Although Eq. (2-23) is general and it is very straightforward to calculate it, the calculation 
is absurdly tedious because J is not a time-ordered operator but a product of ones, and has 
particle-number-changing interactions. In order to renormalize the theory perturbatively, 
we need to calculate the counterterms but it is not manifest what sort of the divergences 
emerge in the effective Hamiltonian before doing actual calculations. It is not obvious either 
whether it is renormalizable or not even in the renormalizable theory. 
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§3. Construction and renormalization of the effective Hamiltonian 

In this section, we give the construction method of the FSTO's effective Hamiltonian from 
Feynman diagrams and discuss the renormahzation of it. First, we introduce the auxihary 
operators F and G for convenience, and then show that the effective Hamiltonian is written 
in terms of them. G can be constructed from the same Feynman diagrams as those for S- 
matrix elements in the covariant perturbation theory. We give the rules for the construction 
of it and show that it is renormalized by the usual renormalization procedure. Lastly, we 
show that the effective Hamiltonian has the divergent terms even if G is renormalized. We 
discuss their role in diagonalizing the effective Hamiltonian. 

3.1. Definitions of F and G 
Let us define an operator 

F = r]H'Jr], (3-1) 

motivated by the first few terms in Eq. (2-23). We will show that the effective Hamiltonian 
can be rewritten in terms of Hq and F. By using Eqs. (|2-19D and (2-20), we write J as 



J=(° ° l + l \ ° I = (1 - r/) + f/(0)r/(f/(0))-V (3-2) 

Therefore, F becomes 

F = r]H'U{0)r]{U{0))-^r] 

= Gv{Um-\ (3-3) 

where we define another operator G as 

G = 7]H'U{0)r] 

= r]T |i/'(0) exp ~i dt H'{t) | r/ 

+ E (-0" / dtj dh--- dt^ H'{0)H'{t,)H'{t2) ■ ■ ■ H'itn) V- (3-4) 



n=l 



As we will shortly see, to introduce G is crucial for using the the diagrammatic rules. An 
advantage for using G is that it is diagonal in the particle number space so that we can 
normal-order the products of G"s more easily than those of J's and H'. Note that G is 
similar to the S-matrix operator but there is one important difference: the upper limit of 
the time integration is 0. 
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We can find that T]{U{t))ri is written in terms of G, by integrating t' from — oo to t after 
sandwiching G between e*^"*' and e~^^°^' : 

-I /"* di! Gii!) 

= -i f dt' r]\H\t') 

J —oo I 

oo 

+ EH1 

n=l 

= v\-i f dt' H'it') 

+ E i-^r^' / dt' / < I dt',--- ""^ < H'{t')H'{t[)H'{t',) - - - H'it'J ^ 

J —oo J —oo J —oo J —oo \ 



^"+^ ' dti / dt2 
-oo J —oo 



dtn H'{t')H'{ti + t')H'{t2 + t') • • ■ H'(tn + t') h 



?i=l 

= v{u{t))v-v, 

where 

By substituting (f/(0)) in Eq. (3-5) to Eq. (3-3), we can write F as 



(3-5) 
(3-6) 







T]-t dt G{t) 

J —oo 

oo r „0 



-1 



F = G 

= |r? + E dt G{t) I . (3-7) 

It is important that F depends only on G. Here, we expand F and G in the order of H': 

oo oo 

-p" = E ^ ^ E ('^■^) 

It is convenient to define 



m=l 



m=l 



(3-9) 



m=l 



m=l 



Also, we solve the differential equation 



dt 



-v{H\t) + H'it)Mt)}ri + ri{H'it) + Mt)"^ H' {t)}r] 
-J'2it) + J^2it)'f , (3-10) 



which is obtained from Eq. ( p-14D . It is easily found that 

jjji=z r dt{j^2{t)-j'i{t)}, 

J —oo 



(3.11) 



under the initial condition Ji(— oo) = 0. From Eqs. (|3-1| ), ( |3-9|) and (|3-11|) , the term 
r]{[JiJi, H'Ji] + h.c.)?7 in Hes reads 

[jUi,vH'Ji]+h.c. 
= [jUi, vH'v + vH'Ji] + [vH't] + jIh't], jUi] 



I 



f dt[T2{t)-Tl{t),T2]+\i.c. (3-12) 



Now we are ready to rewrite the effective Hamiltonian in terms of Hq and F. We finally 
obtain the effective Hamiltonian in the useful form from Eqs. (2-23), ( p-l|) , ( p-9|) and (3-12): 



i^eff = \v [Ho + ^4 + ^ f_Jt [^2(t) - 4W,-^2]} V + h.c. + (3-13) 

All interactions are written in terms of Qn through JF„ to the fourth order in H' . To sum- 
marize, the effective Hamiltonian can be easily constructed once Qn is obtained. 

3.2. Rules for the construction of G 

In this subsection, we give the rules for the construction of G from Feynman diagrams. 
We do not need to use the old fashioned perturbation theory. The knowledge in the covariant 
perturbation theory helps us to find them. 

From Eq. (3-5), we immediately find that 

/oo 
dt G{t), (3-14) 
-oo 

where S is the familiar S-matrix operator. Comparing order by order, one can relate the nth 
order term of 5* to (?„ as 

/oo 
dt Gn{t). (3-15) 
-oo 

It is important to make the difference between rjSnri and Gn clearer. It is useful to rewrite 
them in terms of T-product and compare them in order to find the correspondence with 
Feynman diagrams: 

( POO roo roo 

r]Sr,r] = ^^ / dtj dt2--- dt^ r]T{H'{t,)H'{t2) ■ ■ ■ H'{t^))r], (3-16) 

Til J —oo J —oo J — oo 

Gn= ' r dt, r dt2---C dt^_,r,T{H'{0)H'it,)H'{t2)---H'it^_,))r]. (3-17) 

(n— lj!j-00 J-OO J-OO 

Wick's theorem tells us that T-products of H'{t)'s can be written as sums of normal-ordered 
products of creation and annihilation operators with the coefficients being amplitudes, which 
are given by the same Feynman diagrams. The difference between those arises from the time 
integrations. It is apparent that the role of dt is to give the energy conserving delta 
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function 2Tx5{Eout — Ein) to each vertex in the Feynman diagrams, where Eout is an outgoing 
energy from the vertex and Ein is an incoming energy to it. Because the time integrations 
end at in Gn, the energy denominators {—i){Eout — E^ — ie)~^ appear instead of the delta 
functions. The sign of an infinitesimal constant e must be taken positive to ensure the 
convergence of the integrations at t = — oo. n H'(t)^s become equivalent and cancel the 
factor l/n\ in Sn due to n time integrations, while there are (n — 1) time integrations and 
the factor — 1)! is canceled in G„. But there is a H'{0) which is not integrated, therefore 
there are n terms for each vertex which have {n — 1) products of the energy denominators. 

Let us consider the LF Yukawa modelB as a concrete example. In this model, the inter- 
action part of the Hamiltonian has the form of 



The effect of H'^'^^ which is the instantaneous interaction is absorbed into the fermion prop- 
agator, and we hereafter omit it because we assume that LF diagrams are equivalent to the 
covariant one.0 We define i] as the projection operator for the two-body state which consists 
of a fermion and an anti-fermion. Because the power of g is equal to the number of the scalar 
fields and t] : (odd 0's) : 77 = is satisfied, a product of H'^s satisfy 



for odd m, where the number in the brackets in the superscript means the order of g and 
= {x^ + x^)/\/2 is LF time. JF2 and JF4 become 



H' = gH'^''> + g^H'^'\ 



(3-18) 




(3-19) 




(3-20) 




(3-21) 



Gn is written as 



Gn = T.[ I MnVbliPi)bxip2)vi27r)'S%p,-p2) 

^ J Pi JP2 

-nil MnVd\{pi)dx{p2)Tl{'^7^f5\p,-p2) 
\ Jpi JP2 




xvbUpi)dl{h)d^,{l2)bx2{p2)v{^rff6Hpi + I1-P2- I2), (3-22) 



* Lagrangian and Hamiltonian in this model are given in Appendix 

* It is explained at the last paragraph in Appendix Kl 
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where we define 

6^{p -q)= 5{p+ - q+)6\p^ - q^). (3-24) 

We hereafter call p~ "energy" and {p'^,p±) "three- momentum" . The creation and annihila- 
tion operators of the fermion, b\{p) and bx{p), and ones of the anti-fermion, d\{p) and dx{p), 
satisfy the anti-commutation relations 

{bxM,bUp2)} = {dxM,dUp^)} = 2pt{2nf6\pi-p2)6x,x„ (3-25) 

and the other anti-commutators vanish. The minus signs which are attached to the second 
and third term of the right-hand side of Eq. (3-22) come from the normal-ordered products 
of the anti-fermion operators, but are just the convention. and Ain correspond to the 
amplitude for the fermionic and anti-fermionic one-particle state respectively, and we have 
allowed for the conservation of the helicity. As we discuss in Sec. 3.3.1, these vanish under 
the physical (on-shell) renormalization condition. 

Vn is a potential between a fermion and an anti-fermion in G„ and is constructed by the 
nth order Feynman diagrams which have an incoming fermion and anti-fermion and outgoing 
ones. Since pair creations and annihilations of particles from the vacuum are forbidden due 
to the merit of the LF field theory, we must not include the disconnected diagrams. 

There are three differences between our construction rules of Vn and the usual Feynman 
rules for the general amplitudes: 

• Assign a factor 

- 915 - _~ - _ . , (3-26) 

Pout Pin 

to each vertex, where p'^t (Pin) is the outgoing (incoming) four-momentum from (to) 
the vertex. The energy denominator corresponds to the energy conserving delta func- 
tion in the Feynman rules. 

• Assign an independent energy p~ for each propagator and integrate it. We do not 
impose the conservation of the energies. 

• Attach a total energy difference —{Pi + — P2 ^ ^2) as an overall factor. 

It is obvious that the first two differences come from the domain of the time integrations. 
The origin of the last one is that there are n terms which have different products of the 
energy denominators as mentioned above. For example, in the case of the diagram Fig. 1(a) 
we consider the energy denominators obtained from Eq. ( p-17|) with factors of i: 

—i \ / —2 \ / —i \ 



U -qi -ie - k2 - I2 - - ^2 - P2 - , 
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+ 
+ 
+ 



—I 



—I 



' —i \ f —i \ f —i 



yPi +Qi - ki - ie 



Ji +^2 -Qi -i^J \ki -92 -P2 



le 



(3-27) 



vPi + ?i ~ ^1 —i^J y-i +^2 ~ Q.\ —i^j \?2 ~ ^2 — I2 

The factor {—i)^ comes from three time integrations (— i dx^), while the factor i'^ comes 
from four vertices {—igip'-y^cjyil)) . In general, the total factor of i is where V is the 

number of vertices in the diagram. Since the sum is the combination of removing one from 
four denominators, we can make it a product of all the energy denominators: 



-(Pi +^1 -P2 -k ) 



—I 



X 



vPi +?i -ki - ie) \li +k2 

—i \ / —i 



Qi 



le , 



02 -^2 -h -i^J \ki -q2 -P2 



(3-28) 



This factor corresponds to the ehmination of the total energy conserving delta function. 
Finally, V„ becomes 



-{pi + h -P2 -k) Yl j_ 

diagrams 



~ dk^ f°° dk2 



00 27r 
—i 



27r 



dqi 
27r 



00 



dq2 



'2 

00 27r 



X 



P\ -ki -qi - ie) - + ^2 

X ( i^An)^ 



le 



(3-29) 



where An is the usual amplitude constructed from the Feynman diagrams but each energy 
of the propagators is an independent variable which is integrated outside of An- Therefore, 
An is constructed from the propagators Spik) and Apiq) with each independent energy, and 
the vertex —^'75, and the external fermion line ^(pi, Ai), u{p2, X2), v{li,ai), v{l2,cr2), and an 
integration of four-momentum for each loop. The three-momenta are conserved but energies 
are not. 

Although each propagator includes both the loop and independent energy, it is easy to 
split them. For example, in the box diagram Fig. 1(a) the corresponding energy denominators 
are 



—I 



—I 



—I 



iPi + ?i ~ ki — ie ) ' \li + k2 — qi — ie ) ' \q2 — ^2 ~ ^2 



and 



,^1 -Q2 -P2 -^e^ 



(3-30) 
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If we consider q2 as the loop momentum, by shifting momentum as 



ki^ ki + q2, qi ^ qi + q2 and k2 ^ k2 + q2, (3'31) 

the denominators become 

/ —i \ / —i \ f —i \ 



\Pi + 9i ~ ^1 — ie/ ' + k2 — qi — ie J ^ \~^2 ~ ^2 — ^e/ ' 

whose energies are assigned in Fig. 1(b). Because the loop momenta are not restricted by the 
usual conservation law of four-momenta in the covariant perturbation theory, it is apparent 
that the energy denominators do not depend on them. 

Pi ki P2 Pi ki + q2 P2 



I ^ , ^ ^ I 

I I 



I 



51 f A 92 gi + 92 1 A 92 



♦ i92 9i+92f i 

I I 1 ' 



<i I2 li h 

(a) (b) 

Fig. 1. An example of the one-loop box diagrams. In (a), energies are assigned to each propagator 
individually. In (b), energies are assigned to each propagator but the loop four- momentum q2 
is specified. All the three-momenta are conserved in both case. A continuous line corresponds 
to a fermion one, and a broken line correspond to a scalar one. 



3.3. Renormalization of G 

As mentioned in Sec. 3.2, An corresponds to the usual amplitude. What is different from 
the amplitude is that energies of the propagators in An are independent each other. 

Even if the energies are not conserved, we can renormalize An by the usual prescription in 
the covariant perturbation theory because UV divergences from the integrations of the loop 
momenta are local. They emerge as the coefficients of the polynomial of the other momenta. 
If the zeroth order term in those expansions includes a divergence, it is renormalized by 
shifting masses or a couphng constant. The other divergences are logarithmic and depend 
on momenta, then they must be renormalized by the field renormalization. 

After renormalizing An, we must consider the energy integrations. First, we show that 
the new divergence arises when we integrate the energies of the renormalized An with the 
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energy denominators. In Hamiltonian formalism, it is not so trivial where the field renor- 
malization constants come from because the energies are not conserved. We discuss the field 
renormalization and show that an ambiguity which may cancel the above divergence comes 
from the energy integrations. 

3.3.1. Divergences in G 

Even if An is renormalized, new divergences may 

V„. 

This problem does not occur in the one-particle 
states which correspond to the diagrams like in 
Fig. 2. From the on-shell condition 

pl=pl=m^, (3-33) 

and the three-momentum conservation law 

pt = Pt^ Pi± = m, (3-34) 
for the fermion line, the one-particle states conserve the energies: 

2pJ 2pJ 

The same conservation law is applied to the anti-fermion line. It is obvious from Eq. (|3-14|) 
that G multiplied by the total energy conserving delta function is S. Therefore the one- 
particle states of G are equivalent to those of S. If we renormalize the self-energy part under 
the physical renormalization condition, Ain and A4n in Eq. (3-22) vanish. An example of 
the order g"^ will be demonstrated in Appendix p. 

In multi-particle states, even though the external lines satisfy the on-shell condition and 
the three-momentum conservation law, one can not say that the total energy is conserved. 
We will show that such divergences disappear even in the multi-particle states except for the 
case that the outgoing external fermion line has a self-energy part. 

First, we consider the one-particle irreducible part r{k[, /c2, ■ ■ ■) of the renormalized An, 
where k[,k2,--- are the four- momenta of the propagators except for the loop momentum, and 
satisfy the three-momentum conservation law. Since r{k[, ^2, ■ ■ ■) is covariant, the analysis 
for large four-momentum k'/^'s, which is used in the operator product expansion, is valid 
when we investigate the asymptotic behavior for large /c^~'s.li^ Therefore, after considering 
that the loop momenta associated to k'l^'s are as large as k'l^'s, we regard k'~'s as larger than 



arise from the energy integrations in 



Pi 




P2 




Fig. 2. A diagram of the one-particle 
state. The blob is a connected dia- 
gram. 
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fe-^'s and A;-j_'s. By the power counting, the leading contribution comes from the self-energy 
part. For the scalar self-energy part, the leading behavior of r{k[, ■ ■ ■) goes 

ri^k^, ^2; ' ' ')scalar self— energy ~ ^ ] ^ij^i'^j 

^^Cij{k'+k'r + k'rk'+), (3-36) 

and for the fermion self-energy, 

r{k^^ ^2; ' ' ') fermion self— energy ^ ^ j ^i^ji 

i 

^^c,k'-^+, (3-37) 

i 

where Cij and q may depend on logarithmic factor of /c^'s. It is important that both are 
proportional to k'~^s at most. If we multiply r{k[,k2, ■ • •) by the energy denominators and 
integrate over the energies, it seems that it behaves as logarithmically divergent. But such 
term must not exist when the energies are conserved, that is, such divergence is not possible 
in S, so it must be proportional to the total energy difference. For the scalar self-energy in the 
internal line, we regard the blob in Fig. 3(a) as r{k[, k2, ■ • ^scalar self-energy- The contribution 
from it is 

dk^ dk^ -i 
J-oo 27T 7-00 2n "' K - k'^ - ie" ^ ' ' '^''^''^ -'f— -gy 

~ — = -{^^(^1, 92) + iq^{q2 - Qi) log(g'+yl") 

Q2 ~ Qi ~ 

+{q2 — (finite terms)}, (3-38) 

where 77(gi,g2) becomes the usual scalar self-energy if we replace the energy denominator 
with the energy delta function, and /l^ is the UV cutoff of the energy. The second and third 
term are proportional to {q2 — qi ) by the above reason. g+ and g'+ are proper longitudinal 
momenta. The divergence only arises in the second term because Eq. (3-38) is regarded 
as Taylor expansion in (g^ — qi)- The first term is the ordinary renormalized finite term. 
Although the second term is divergent, it is found that it behaves as logarithmically divergent 
at most by Lorentz covariance and the power counting . The third term converges because 
the second term is logarithmically divergent and a differentiation by the energy decreases 
the power of the divergence by one. 

Although Eq. (3-38) includes the divergence, it vanishes when it is inserted in the internal 
scalar line. Let us consider scalar propagators and energy denominators in vertices around 
the blob in Fig. 3(a). It is crucial that the second term in Eq. (3-38) does not depend q^ 
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and ^2 • The total contribution from the second term is 

dqi f°° dq2 -i a f \ 

^F{q2) 



3 27r J-oo 27r — — q2 — ie 
X I -iq'^{q2 - qi) log(g'+yl-) \ Apiq,) — 



. ^2 - 9i - «e J q^ + k2 - fci -le 

= AF{k, - ks)ei-qt)q+ \ogiq'+A-)AF{k, - k2)d{qt). (3-39) 

But the three-momenta are conserved [q^ = q2,qi± = q2±), so it is obvious that Eq. (3-39) 
vanishes. 

The fermion self-energy is logarithmically divergent like Eq. (3-38). For the internal 
fermion line in Fig. 3(b), the contribution from the divergent part is 

dkn dkn 



^3 



SFih) 



2% J-oo 2ti k^^ + q2 — — ie 

where k^ is a proper longitudinal momentum. Unlike Ap^q), 5'_p(A;) has a term which does 
not depend on the energy but is proportional to 7^: 

.-^^;^^2l±p:_Lii^j±±^. ,341) 

2/c+ k^ — + te 

Fortunately, 7+ has the property that (7''')^ — 0, so the first term has no effect in Eq. (3-40). 
The second term is discussed in the same way as in the scalar case, and then Eq. (3-40) 
vanishes. 

Note that the above mechanism only works for the internal lines. The diagrams which 
have the self-energy in the external fermion line may have the same divergence. We consider 
the incoming external fermion line like in Fig. 3(c). The contribution from the divergent 
part becomes 

iy-t q,^k,~-k,-J ^^'^^ {^r^i3-^7-(^r -Mk^A-)}uiP2,x2) 

= SF{qi + k2)e{-pt) \og{k+A-)u{p2, A2), (3-42) 

where all the three-momenta are conserved. Since the incoming line has the positive longi- 
tudinal momentum pt > 0, this vanishes. In the case of the outgoing external fermion line 
like in Fig. 3(d), the contribution from the divergent part becomes 



00 f j ~) j 

' u{pu Ai) - . i7+(pr - ^r) logik^A-) SF{k,] 



27r [pi — ki — ie J ki — k2 — qi — ie 

= u(pi, Xi)SF(k2 + qi)9(pt) \og(k+A-). (3-43) 

Since the outgoing hne satisfies p^" > 0, this does not vanish and is logarithmically divergent. 
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(c) (d) 
Fig. 3. (a) is a diagram which includes a scalar self-energy in an internal scalar line, (b), (c), 
and (d) are ones including a fermion self-energy part in an internal, incoming-external and 
outgoing-external fermion line, respectively. In all the case, the three-momenta are conserved. 
The blobs are some one-particle-irreducible graphs. 



3.3.2. Field renorm.alization 

The counterterms for the fields have the same behaviors as those in Eqs. (3-36) and 
(3-37). The energy integrations of those yield not logarithmic divergences but ambiguous 

constants. 

First, let us consider the case of the scalar field. We define the field renormalization 
constant as Z^. The corresponding counterterm in our rules is 

^ - ,! ■ (1 - Z,)A-/{k)- (3-44) 

-OO ZTT g'2 — A; — ^e k — qi — le 

where the three-momenta are conserved [k^ = qt = qt, = qi± = q2±) and we omit the 
scalar propagators attached on both sides of it. This is the usual field renormalization and 
cancels the logarithmic divergence which comes from the integration of the loop momentum. 
Carrying out the energy integration, we obtain 



cqt + — = -^f\<1i) 

?2 



(3-45) 



where c is an ambiguous constant which depends on the way of taking the limit k~ — > ±oo. 
More generally, the contribution from n counterterms connected by the scalar propagators 
is 



Fig. 4(b) 



dki dk2 dk^ /■°° dq2 dq^ dq, 



27r J-oo 27r J-oo 27r J-oo 27r J-co 27r J-oo 27r 
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Qn+1 



k„ — ie 



—I 



^2 Q2 



le 



-AF{q2 



92 -h 



It 



—I 



Qn+1 -Qi - «e 

where all the three-momenta are conserved. is defined as 

Sn(g+) = a„^(g+) + 6„0(-g+), 



ki — Qi — ie 
(3-46) 



(3-47) 



where a„ and bn are ambiguous constants because there is an ambiguity of the order of the 
energy integrations.0 If we insert it in an internal scalar line, the contribution becomes 



Fig. 4(c) 



^ d(h_ 



27T ^3 Qn+l 



Apiqn 



+1> 



Qn+1 Ql 



It 



^^'^^^ Qi + k2 — ki — le 
Apik, - k^)e{-q+^^){l - Z^rSn{qt)qtAF{ki - k2)9{qt) 



+ 



dq 



00 2tt k^ — k^ 



q^ — It 



;i-Z3)M^(g)- 



q' 



k. 



dq~ 



—I 



00 277 /C4 — 



-{l-Z,rAF{q) — 



2 ^1 

i 



. (3-48) 



q —le' ' "q +k2 —ki —le 

This does not depend on Sn{q~^) and is equivalent to a free propagator multiplied by the 
constant (1 — The mechanism that vanishes is the same as one that the loga- 

rithmic divergence in the self-energy does. We take the sum of the diagrams like in Fig. 4(d) 
and obtain the result 

dq~ —i 



i + Y^{i-z,r 

V, n=l y 

1 r°° dq~ 



-00 27r — 



q' 



—I 



Z-i J -00 2% ki — fco 



q — le 



-Ap{q) 



- le 
—i 



-AM 



—I 



q + k2 — ki — ie 



q + k2 — ki — ie' 



(3-49) 



where the three-momenta are conserved. This property is the same as that of the usual field 
renormalization. 

Similarly, we can discuss the fermion field. If we connect n counterterms for the fermion 
field with fermion propagators, the contribution becomes 



Fig. 5(a) 



dqi 



dq2 



2tx J-oo 2tt 



dqn 

27r 



CXD 



dkr, 



'2_ 

■00 2tt 



00 



dkr. 



■3_ 

00 2tt 



27r 



See Appendix 
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(a) 



(b) 






92 qi 




(c) 



+ 



(d) 

Fig. 4. Diagrams which contribute to the renormahzation of the scalar field. X means a coun- 
terterm for the scalar field, (a) is a diagram which have a counterterm for the scalar field 
between scalar lines, qi and q2 are four-momenta assigned to the external scalar lines and k is 
one for the counterterm. (b) is a diagram which have n counterterms connected by the scalar 
propagators. In (c), (b) is inserted between two vertices, (d) is an infinite sum of the diagrams 
(c). In all the case, the three-momenta are conserved. 



1+1 Qn 



-{I - Z2)S-p\q^)- 



{l-Z2)S-p\q2) 



?2 ~ ^2 



-SF{k2 



-i 



le 



k2 -Qi -ie 



= (1 - z^y 



fn{kt)l^ + 



-S-F\k^] 



Qi — ki — ie 
(3-50) 



kn+i ki ie 

where all the three-momenta are conserved and Z2 is the field renormalization constant for 
the fermion field. The function fn{k^) is defined as 



fn{k+)=cJ{k+)+dJ{-k+), 



(3-51) 



where c„ and d„ are ambiguous constants due to a similar ambiguity. 

When we insert Fig. 5(a) into the internal fermion line, the contribution from it is given 

by 



Fig. 5(b) 



°° dki 



271 



dk. 



n+l 



—I 



27r q2 +P2 - kn+i - ^t 



Spikn+i) 
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X(l-Z2)' 

x-S^(A;i)- 



fn{kt)r + 



— I 



ki -Qi -Pi -ie 



+ 



oo 2% q2 + P2 — k — ie 

°° dk- -i 



k — Qi — Pi — ie 



1 - Z2TSFik)- 



-, (3-52) 



00 



271 ^2 + P2 — k~ — ie k~ — Qi — Pi — ie' 

which is a free propagator multiphed by the constant {I — Z2)"' and the energy denominators. 
An infinite sum of the diagrams becomes 

^ r ^ =^ SAkh ^ . (3-53) 

Z2 J -00 27r q2 + P2 — k — le k — Qi — Pi — le 

As a result, fn{k'^) does not affect the fermion propagator. 

For the incoming external fermion line, the contribution from the n connected countert- 
erms is 



Fig. 5(c) - (1 - Z2nfn{kt)9{-kt) + A) 



—I 



Qi +Pi - ki - ie 



[l-Z2ru{kiA) — 



—I 



(3-54) 



^1 +Pi —ki — ie' 

where kl — and ki > 0. Of course, all the three-momenta are conserved. In this case, 
fn{k~^) does not affect it. The infinite sum becomes 



■^u{k, — - \_ -, 

Z2 Qi + Pi - k - ze 



(3-55) 



which is I/Z2 times the tree external line, k'^ is the momentum of the external fermion. 
For the outgoing external fermion line, the contribution is 



Fig. 5(d) = (1 - ^2)"[/n(^++l)^(^:+l) + l]u{kn+U A)-3 z 

kn+1 -Qi - Pi - «e 



{l-Z2r{Cn + l)u{kn+l,X) 



—I 



kn+1 ~ Ql ~ Pi ~ 



(3-56) 



where /c^^^ = and k^+i > 0. The infinite sum becomes 



l + ^(l-Z2)"(c„ + l)U(fc,A) 



n=l 



k — Ql — Pi — ie' 



(3-57) 



which is the tree external line multiplied by the constant. Note that c„ may not be equal to 
zero. 



20 



«2 fc2 fci '^^X kn+i ■■■ fc2 A:i 

P2 / \ Pi 



rK *-{ 



(a) (b) 

\ I 

9l V , . , , ; , , V gi 

\ Kri+1 K2 Kl fCn+i ^2 / 

r K K ^ K ^ external external ^ K K ^ K ^ ( 

(c) (d) 
Fig. 5. Diagrams which contribute to the renormahzation of the fermion field. X means a coun- 
terterm for the fermion field, (a) is a diagram which have n counterterms connected by the 
fermion propagators. In (b), (a) is inserted between two vertices. In (c), (a) is attached to the 
vertex and connected to the incoming line, and k\ = m? and /c^ > is satisfied. In (d), (a) is 
attached to the vertex and connected to the outgoing line, and kn+ii = m? and A;^_^^ > is 
satisfied. In all the case, the three-momenta are conserved. 

3.3.3. Summary of the renormalization of G 

We have shown that UV divergences from the loop integrations are renormalized by the 
usual procedure in the covariant perturbation theory. After renormalizing An, we carry out 
the energy integrations. As a result, the new UV divergence and ambiguous constant c„ 
arise only from the diagrams that the fermion self-energy and its counterterm are inserted 
into the outgoing external line. 

We fix Cn so that it can remove the above divergence. This is always possible because 
both live in the same place and are closely related. Although it is not the only way to fix 
Cn, we consider that it is natural for both to cancel each other. 



3.4. Divergences in H^s 

In Sec. 3.3, we have shown that G is not divergent after renormalizing it by the usual 
prescription in the usual perturbation theory. In this subsection, we show that i^efr is 
divergent even if G is finite. 

From Eqs. ( ^-131 ), ( p-20| ), and ( p-21| ), the effective Hamiltonian is written as 



+ 7/° ^^a;+[6;2(x+)-d(x+),6;2])r/ + h.c. + 0(/), (3-58) 
4 J —00 J 
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to the order of g^. We can renormalize Q2 and ^4, but new divergences arise from the product 
of Q2- One of the examples is a product of two one-scalar-exchange diagrams in Fig. 6. The 
corresponding term in Q2 is 

GT^-iJ2J2J2^ / / u{pi,\i){-g-fr,)u{p2,\2)v{l2,(J2){-gib)v{lu(Ji) 



Ai A2 ""1 <^2 



X [A^ik - h)e{lt - It) + AAlh P2)0{pt - pt)] 

xr)bl{pi)dl{h)dMK{p2H'^7rfS%Pi + h-p2- I2), (3-59) 

where we have integrated the energy of Aplq). The product of two ^|^'s included in the 
third term of Eq. (3-58) is 



J— 00 



log ''i E E i I i i v/^a-Mi-') 



A cr 



^^j e{x-x') /\og§ ^ i\ ^ e{x'-x) / log^ ^ 1 



1 — x' \1 — X X I x' \ X 1 — X 



xr)bl{pi)dl{h)d^{l2)bx{p2H'^7TfS\pi + I1-P2- I2) + (finite terms), (3-60) 



where Ai_ is the cutoff of the transverse momen- 
tum of the external fermion. x and x' are defined 
as 

Pi , P2 



Pi 

— r 



a; = — r — TT,x = — r — -r, (3-d1) 1 

Pl+lt P2+12 ^ ^ — 

respectively. Eq. (3-60) is logarithmically divergent h h 

as yl_L ^ 00. A Feynman box diagram, which is Fig- 6. The one-scakr-cxchange dia- 

finite by the power counting, consists of the sum of ^^^^^ contributes to ^2- 
the various time-ordered diagrams which may be divergent individually. The above product 
is one of such and is logarithmically divergent. 

It is important to recognize that H^s should be divergent. The reason is understood by 
the perturbative diagonalization of it: 

^pert = \H [Ho + Qa 



2 

+ 



^6^2/° dx+g2{x+) + r dx+g2{x+)g2 + h.c.]\a) + 0{g'), (3-62) 
z J —00 Z J— 00 ) 



where \a) is an eigenstate of Hq. Because the logarithmic divergence in Eq. (3-60) comes 
from the integral of the three-momentum of the intermediate state, the divergences of 
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I J-oo dx^Q2{x^)G2 and of /_oo '^^^^2(3^^) have the opposite sign, and they cancel each 
other. In general, using the eigenstates of Hq and its complete set, we obtain the matrix 
element 



Z J —00 L J —00 

E(«l^2|7)(7l^2|/3)^ (-—^—— + 



7 



2 -pp-ie p-^ 



= E(«l^2|7)(7l^2|/?)^^^ J^'-.J^^ „ y (3-63) 

in which the power of the intermediate energy p~ is decreased by one, and thus the power 
of the transverse momentum of the intermediate states by two. 

It is important that although these divergences must be canceled by adding the artificial 
counterterms in the other similarity methods, they will automatically arise in the higher order 
terms in our method. As shown in Appendix 0, if we diagonalize the effective Hamiltonian 
to the order of g^, the eigenvalue is divergent. In general, such divergences are mainly related 
to the box diagrams and arise not only in this case but also in the case of TD approximation 
in (3 + 1) dimension.© Eq. (3-63) says that such divergences are canceled if we include 
order terms, then the divergences in the eigenvalue are weakened. However, new divergences 
arise in diagonalizing the effective Hamiltonian due to new interactions. We expect that 
they are canceled by the higher order interactions because the exact eigenvalue should not 
depend on the cutoff and the similarity transformation does not change the eigenvalue. 



§4. Summary and discussions 

In this paper, we have shown that the effective Hamiltonian, which is obtained by the 
FSTO's similarity transformation in the particle number space, can be written in terms only 
of G (or F) to the fourth order in H'. 

To introduce G is crucial for constructing the effective Hamiltonian more easily than in 
the traditional way especially in the higher orders. G has a favorable property that it is 
diagonal in the particle number space. Since it also has the same form of the formula of the 
S-matrix operator, we can use the Feynman diagrams and the rules mentioned in Sec. 3.2 in 
the LP Yukawa model. By using the knowledge in the covariant perturbation theory, we can 
avoid the complexity in calculating it and make the renormalization procedure transparent. 
Note that our construction rules are a little bit different from the usual Feynman rules 
because energies are not conserved in each vertices. 
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The divergences due to the integrations of the loop momenta are renormahzed by the 
famihar prescription in the covariant perturbation theory. The divergences due to the energy 
integrations can be canceled by terms which come from the ambiguity of the counterterms 
for the fields. 

The mechanism of the cancellation of the extra divergences from the energy integrations 
is valid only in the LF field theory. It will be applied not only to the LF Yukawa model but 
also to the other LF models. Although we can not use it in the ET field theory, we expect 
that more precise integrations and treatment of the limit are necessary for the cancellation 
of the divergences. 

Although G is constructed by Feynman diagrams and renormahzed, the effective Hamil- 
tonian has divergent terms which are written in terms of products of the renormahzed those. 
It is very important that such terms do not need the counterterms but work as those which 
cancel the divergences in diagonalizing the effective Hamiltonian. Although the exact eigen- 
value should not depend on the cutoff, our method is perturbative, so it is not possible to 
treat nonperturbative divergences in diagonalization and to get the exact eigenvalue. It is 
highly desired to find the nonperturbative counterterms for it. When we find the low energy 
states, our method is useful in small coupling region enough to allow one to ignore the cutoff 
dependence if the cutoff is much larger than lower eigenvalues. 

There is also the problem of vanishing energy denominator. Here, we only consider the 
LF Yukawa model without a massless particle, so that we may avoid this problem by —ie in 
denominators. If there is a massless particle like in QCD, we must be attentive to it. 

We proved here that the effective Hamiltonian can be written in terms of G to the fourth 
order by the explicit calculation. But the general proof for higher orders is lacking. Although 
we do not know how the general proof goes, we think that it is likely that this feature persists 
to all orders. Recently Hansper proposed a nonperturbative approach for the FSTO's method 
in LF field theory. © Although he applied it only to the parton distributions, it is very useful 
to solve the mesonic bound states in QCD if it is applicable to our method. 

In Appendix |B|, when the coupling constant is large enough, the lowest-energy state has 
the positive binding energy in the numerical calculation, but we do not regard it as a bound 
state because it depends on the transverse cutoff yl_L. As mentioned above, our method is 
valid in small coupling and we expect that it is improved in higher order calculations. We 
are now extending the present work to the next order in order to confirm the cancellation of 
the cutoff dependence and to get bound states. 0^ 

It is interesting to note that this method might also be applied to the similarity trans- 
formation in momentum space because Eq. ( ^-13 ) does not depend on the choice of rj. If so, 
to get the effective Hamiltonian would be much easier than the traditional way. 



24 



Acknowledgement 



The author is grateful to K. Harada for helpful discussions and encouragement, and also 
thanks A. Okazaki for his useful work. 



Appendix A 
Light-Front Yukawa model 

In the LF framework, we use = (x'^ + x^)/\/2 as a time variable instead of t. The 
Lagrangian of the this model is 



2 



+ 'ipii'j^df, -m + ig'j5(j))ip. 



(A-l) 



By LF quantization with Dirac's quantization method, dynamical fields must satisfy the 
following commutation relations 



1 



6{x -y )5^{x± -y±), 
6{x' - y')6^{x± - y±)A+, 



where ^(x) is the dynamical part of the fermion field ^/'(x): 

= A+ij{x). 

A± are the projection operators of the fermion field and defined as 

1 



A, 



V2 



7V- 



From Eq. ( |A-1| ), the familiar Legendre transformation gives the LF Hamiltonian 

H = P- = Ho + H', 



where 



Ho 
H' 



dx-d'x^ { ^{d^<PY + lfi'<P' + 



J dx d'^x± 



^2 



^"^7°750 ^^^ {ia±-d± + m(3)^ - h.c. 



+ - 



(0r)(r0) 



(A-2) 
(A-3) 

(A-4) 
(A-5) 



(A-6) 



(A-7) 



(A- 
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and 



In the interaction picture, they become 



1 



—9^ + 



(A-9) 



(A-IO) 



H\x~^) = J dx d'^x± < —ig4''j5(p4' + ^ ' 



+HcTix 



(A-ll) 



where all the field operators are defined in this picture. Note that ip{x) is not the one in 
Eq. ( |A-1D but a new field made from ^(x) in this picture, and the solution of the free Dirac 
equation: 0* 



1 + 



1 



{ia±-d± + m(3) 



ax) 



(A-12) 



Hct{x~^) is the part which corresponds to the usual counterterms obtained by shifting the 
parameters. 

We assume that the LF diagrams are equivalent to the covariant ones. More concretely, 
in the LF framework the fermion propagator 



Sp{k) = / rf^x(0|T+(^(x)V^(0))|0)e 



ik-x 



^ — m + ie 2k~^ 



7 



(A- 13) 



includes the noncovariant term. But in the first two diagrams in Fig. 7, the contributions 
from the second term are naively canceled by the vertex 



1 



- ^9 7 



1 



2{ki + q+) 2{kt-qt 



(A- 14) 



in the last diagram in Fig. 7, which comes from the instantaneous interaction in the Hamil- 
tonian. Therefore, we regard the fermion propagator as the first term in Eq. ( |A-13| ) and 
covariant effectively, and omit the second term in Eq. ( |A-13| ) and the vertex Eq. ( [A-14| ) 
together. As a result, the scalar and fermion propagator are given by 



— IJ? + ie 



(A-15) 



— m + ie 



(A- 16) 
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fc2^ 





Fig. 7. The first two diagrams are examples which include an internal fermion line. The propaga- 
tors include the noncovariant term. The last diagram is the contribution from the instantaneous 
vertex included in the LF Hamiltonian. 



Appendix B 

Explicit calculations and numerical result 

We calculate the ground state energy numerically to the second order in g. Of course, 
the calculations to this order are not so different from the other methods; 

0) 

the advantage 

of the present formulation becomes apparent in the higher orders. S-*. The reason why we 
present the second order calculation here is to demonstrate some features of our method and 
to clarify what would be expected in the next order. 

From Eqs. ( PT3D and (|3-20|) , the LF effective Hamiltonian is 



(B-1) 



to order. Also, we add the flavor of fermions. It is easy to estimate Eq. from 
Feynman diagrams by using our rules. 

Since possible graphs are the same as those needed in constructing the S-matrix, we can 
immediately imagine those which contribute to Eq. ( p-l| ). The Feynman diagrams associated 
with the effective Hamiltonian to the second order are shown in Fig. 8. 

Terms which should be renormalized are only fermion self-energy terms M.2 and in 
this order. Although we do not give the rules for those, they are much the same as those for 
Vn- The self-energy terms correspond to the first and the third graphs in the second line of 
Fig. 8 and are written as 



^2 



27r \pi — k — ie I \k 



P2 -^e. 



(2vr)^ 



+ 



{l-Z2,Y'^SFl{k)-z6mP \ui{p2,X) 



(B-2) 



27 



+ 



+ 



+ 



+ 



T 



>-< 



Fig. 8. Feynman diagrams which are necessary to construct the effective Hamiltonian to the second 
order in g. The first diagram is the free part. The next four diagrams are self-energy parts 
and X mean the counterterms for the fermion mass and field. The last two diagrams are the 
one-scalar-exchange part and the fermion annihilation part. 

where dmf'^^s and (1 — Z2j)*'^'''s are the masses and field renormalization constants of order 
respectively, and the on-shell condition {pf = p\ = mf) for the external fermion lines is 
satisfied. Although there is no energy conserving delta function, pi = P2 is satisfied due to 
Pi = pt^ Pi± = P2± and the on-shell condition. Therefore 

-(pr -P2 



Pi — k — ie J \k — P2 ~ , 
1 1 



Pi — k — ie k — pi — it 

= 2m5{p'{ - k~). (B-3) 

This is what is said in the beginning of Sec. 3.3.1. Integrating we obtain 

<3t' = Y.Y.j A)rf A) vhUp)hx{p)ri. (B-4) 



A 

where 



Sf\p) = ^ J ^,AMi-91n)SF.{p - q)h91,) + Sm^\ (B-5) 

and p"^ = mf. The counterterm for the fermion field in Eq. ( |B-5| ) vanishes by using the 

(2) 

Dirac equation. Obviously, El (p) is the same as fermion self-energy terms in the covariant 
perturbation theory. It is convenient to impose the physical renormalization condition 

u,{p, \)uf\p)u,{p, A) = (p2 = m^), (B.6) 
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in order for the fermion masses to be the physical one. For anti-fermion, the same argument 
is apphed. 

The rest are the scalar exchange part and the fermion annihilation part which are repre- 
sented by the last two graphs in Fig. 8. After integrating the energies, we obtain 

4 



Ai A2 o"! o'a 



1 1 



X 



+ 



Al A2 '^l <^2 -^Pl -^^2 Jh JI2 

XUi{pi, Xi)'j5Viik, Cri)Vi{l2, Cr2)75Mi(p2, A2) 
1 1 



(B-7) 



X 



+ 



(pi+hr-fi^ (P2 +12^-1^' 



X Vbjx, (Pi ) dU ik) di<T2 ik) hiX2 {P2)v, 



(B- 



which are all the second order interactions in the LF effective Hamiltonian. Of course, since 
it does not have the particle-number-changing interactions, the eigenstate is a pure two-body 
state. 

We set the total transverse momentum P± to for simplicity. The eigenstate may be 
written as 



pi POO p2tc 2 



Al A2 



{2-Kf\x{l-x) 



^g^(.n-Ai/2-A2/2V^..(^^ Al, A2, m)4i(Pl)4A2(P2)|0), (B-9) 



7t 



where 



p+ = xP+, P2 = (1 - x)P^, (pi_L - P2±)/2 = 



K 



tan(/7 = ^. 



(B-10) 
(B-11) 



m is the eigenvalue of the third component of the total angular momentum J3. In the LF 
coordinates, the squared total angular momentum is not a good quantum number, but 
the third component J3, the helicity, is. 

We discrctizcd x with L equally- spaced points, in the numerical calculations of diagonal- 
ization. We put a transverse cutoff A± for k and, also used the following variable 



z — 



(B-12) 
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instead of k and discretized z with equally-spaced points because the wavefunctions are 
sharp for /t ~ 0, but fiat for k, ~ A±. 

The results are shown in Fig. 9 and Fig. 10. We set the parameters L = 10, = 30, 
m = 0, the fermion mass mi = 1.0 GeV, the anti-fermion mass m2 = 1.0 GeV and the scalar 
mass fi = 0.01 GeV in all the cases. We define ag = g"^ /'i'K. 

Fig. 9 shows Ug dependence of the binding energy of the ground state for various trans- 
verse cutoff A^. (a) is the case that the fermion flavor is different from the anti-fermion one, 
that is, exchiding the fermion annihilation part, (b) is the case of inchiding it. In (a), the 
binding energy is ahnost independent of ylj_ in < 1.5, but the ground state is not a bound 
state because it is slightly negative. Even though the binding energy becomes positive and 
larger as grows, we can not say that the ground state is bound because it apparently 
depends on for large Ug. This means that the binding energy depends on even if ag 
is small. If we consider the perturbation theory of the effective Hamiltonian, which treats 
^(^2 + ^1)^/2 as the interaction, with small ag, the leading term of the eigenvalue which 
depends on the cutoff corresponds to Eq. (3-60) times (—1). As mentioned in Sec. 3.4, we 
expect that the dependence on A±_ is weakened in the next order calculation, (b) has stronger 
dependence on A±_ than (a). 

Fig. 10 shows yl_L dependence of the binding energy of the ground state for = 3 with the 
fermion annihilation part. The reason for its behavior is not that the LF effective Hamilto- 
nian needs the renormalization of the scalar mass which corresponds to a fermion loop. The 
perturbative analysis shows that the fourth order term of the eigenvalue is logarithmically 
divergent hke in Eq. (3-60). We consider that it is canceled by adding the fourth order terms. 



Appendix C 
Ambiguity in the integration 

In this appendix, we explain that there is in general an ambiguity in the integrations 
like that appearing in Eq. (3-46) by demonstrating an explicit example. Let us consider the 
following integration: 

7= lim lim r ^ r [ Zi . (C-l) 

Ai^ooA2^ooJ-Ai 27r J-A2 2tt ki — Qi — ie Qi — q2 — i€ ^2 ~ ^2 ~ 

If ki and k2 are finite, we can take the integration contours which have a infinite radius in 
the lower half plane for and in the upper one for respectively. By the residue theorem, 
/ becomes 

^ ki +k2 i 1. . d9i 662 1 



ki — k2 —16 yli-^-oo yl2-»oo 



lim lim i f — ^ / — ^ 2 — --; — t-^, (C-2) 

ii-^ooyi2-»oo Jo 27r Jo 27r 1 - 4^e<^2-^i) ^ ^ 
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(a) _ _(b) 
Fig. 9. ag dependence of the binding energy of the ground state for various A±. ag is defined as 
/At:. The calculations are done for L = 10, = 30, m = 0, mi = m2 = 1.0 GeV, and 
fj, = 0.01 GeV. (a) is the case of excluding the fermion annihilation part and (b) is the case of 
including it. 



Ai dependence of the bound state (including annihilation) 




Ai [GeV] 

Fig. 10. A± dependence of the binding energy of the ground state including the fermion annihila- 
tion part. The calculations are done for L = 10, N = 30, ag = 0.716, m = 0, mi = m2 = 1.0 
GeV, and n = 0.01 GeV. 



where we keep the radii finite in the second term because / have various values depending 
on ways to take the limits. For example. 



1. vli — > oo before A2 00, 



2. A2 ^ 00 before Ai — > 00, 



2^' k^-k2-ie^l- ^^'^^ 



^ rv-^ ni'2 
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3. yli = yl2 — >■ OO, 



ki+k^ i ^^^^^ 



2 ki — k2 — ie 8 

The other hmits may yield the other constants. Note that such ambiguity arises from the 
terms whose total dimension of the variables of the integrations is zero. 
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